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Abstract. We apply a method recently devised by one of the authors to obtain 
precise analytical formulas for the spectrum of quartic and sextic anharmonic potential. 
Due to its general features the method can be applied with minimal effort to a wide 
class of quantum potentials thus allowing very promising applications. 
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Recently one of the authors has developed a method which allows to obtain 
analytical approximations to a certain class of integrals with arbitrary precision PQI21- 
It has been proven that the expansion proposed in pQ converges uniformly and that it 
yields excellent results already to first order. Such method has been applied to a large 
class of problems in Classical Mechanics and General Relativity, allowing to obtain 
simple analytical formulas which work even in the non-perturbative regime. 

In this letter we deal with a further application of the method, showing that it 
is possible to obtain accurate analytical expressions for the spectrum of a quantum 
potential by using it in conjunction of the WKB expansion. Although our calculations 
are obtained for the anharmonic potential, the procedure that we suggest applies with 
minor changes to any potential. 

We will use jH] as our main reference: in that paper Bender and collaborators show 
how to solve the one- dimensional eigenvalue problem for analytic potential to all orders 
in the WKB approximation. They obtain a recursion relation for the WKB corrections 
and they are able to sum the WKB series in two special cases, recovering the exact 
energy eigenvalues. 

However, in the general case of a potential V(x) the calculation of the WKB 
correction to a given order can be difficult and, depending on the potential itself, it is 
possible that no analytical expression can be found for the WKB integrals. By applying 
the method of £Q to the anharmonic oscillator to order h 4 we show that one can obtain a 
precise analytical approximation to the entire spectrum of the potential. The precision of 
the approximation depends both on the order to which the WKB expansion is considered 
and on the order to which our method is applied. 

Let us briefly review the method by considering the integral J\ = 



E = V(x±) is the energy. In the spirit of the Linear Delta Expansion jl] we interpolate 
the potential V(x) as 



where Vq(x) is a potential introduced by hand, which depends on one or more arbitrary 
parameters§. V&(x) reduces to the full potential for 5 = 1. We want to perform an 
expansion described by eq. (0) without moving the inversion points; for this reason we 
ask that x± be the inversion points also of the potential Vq(x). As a result, the energy 
E that the particle would possess if it was moving only in the potential Vq{x) will be 
given by: 




x± are the classical turning points of the potential V(x) and 



V s (x) = V (x)+S (V(x)-V (x)) 



(1) 



E = V (x ± ) . 
We therefore write the integral as 



(2) 




(3) 



§ In the following we call A these parameters. 
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where 

E-E -V( X) + V ( X ) 

E - Vo(x) K 
We treat the term proportional to 5 as a perturbation and perform an expansion 
in 5. Such expansion will converge uniformly when |A(x)| < 1 in the region x_ < x < 
x + P]. This condition selects a particular region in the parameter space, i.e. A; however, 
maximal convergence is achieved when the Principle of Minimal Sensitivity (PMS) [3] 
is used, i.e. when the condition dj\jd\ = is enforced. Notice that if the potential 
Vo(x) is chosen appropriately it will be possible to calculate analytically each term in 
the expansion. 

We now come explicitly to the WKB integrals. We perform our analysis to order 
h 4 and use the recursion formulas of ■3j to calculate the spectrum of the anharmonic 
oscillator. Following [Oj we write the integral in a fashion where the singularities in the 
integrand are integrable. The price that we pay is that we have to introduce derivatives 
with respect to the energy. 

We therefore write the WKB condition as 

A(E) = (n + 1/2) (5) 

V 2m 

where A(E) to order 0(h 6 ) is given by 



We have defined: 



2-' + 



J^E)= / JE-V(x)dx (7) 
J 2 (E)= ( X+ V " {X) dx (8) 



E - V(x) 



, .- 7 V"(xf - 5 V'(x) V" ... , , 
Js{E)= / 1 ; , 1 ; y -^dx (9) 



r "'(x) 



'E - V(x) 

where x± are the classical turning points. The spectrum of the potential V(x) can be 
obtained by solving Eq. (JHJ). 

So far we have used general formulas; let us consider now the specific case of 

2 2 .4 

the anharmonic oscillator V(x) = x — |- We choose the interpolating potential 
Vq(x) = \ {muj 2 + A 2 ) x 2 . The method of [T] allows to calculate the integrals in Eq. (JBJ) 
quite simply. Indeed to order 5 10 and using the value of A obtained by the PMS to first 
order one obtains: 

3/2 3 10 

MO = 7, t* E ^ c (io) 

(5 + 8 C)~ «=o 
= ^^(3 + 2C)» I cf 

(21 + 36 C + 16 C 2 )~ n=o 
= m »^(99 + 48 < + 56<»)i - , , 

\/2 (2 (363 + 564 C + 360 C 2 + 224 C 3 )~ «=o 
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Figure 1. Error over the integral Jl for m = u> = % = = 1 as a function of the 
energy. 



where we have introduced the variable ( = — ^ 2 and x± = ±A are the classical turning 
points. The coefficients c$ are pure numbers, which we will not display here. 



{exact) 



x 100 as a function of the 



In Fig. n w e plot the error defined as 5, 
energy. We exploit the relation between the amplitude and the energy, i.e. E = V(A), 
to express ( as a function of the energy: 



c 



2 4 

m u 



4/i-Bn 



l + \ 1 + 



AfxEn 



(13) 



Ji 



[exact) 



are the exact integrals. We have assumed the values m = tu = h = n=l. 
We observe that the error is smaller at lower energies and flattens around values of 
the order of 10~ 5 at large energies. Such result can be further improved by applying our 
method to higher orders in 5: indeed it was proved in pQ that the method is convergent 
and therefore the precision of the calculation can be easily improved. We stress that 
the behavior of this error is somewhat complementary to the one expected for the WKB 
expansion, which is known to work better for the highly excited states. This is of course 
a nice feature, which allows to obtain a higher precision in our calculations. 

Once Eqs. (jTOJ) , (fTTj) and (|T2|) are substituted back in Eq. (JHJ) one obtains a rather 
complicated expression, which still needs to be inverted in order to give the energy as a 
function of the quantum number n. However this task is easily accomplished by Taylor 
expanding Eq. (jUJ) around ( = (which corresponds to the higher part of the spectrum) 
and by then using eq. (fl~3|) . 

We obtain the approximate analytical formula for the spectrum of the quartic 
anharmonic oscillator: 



E as 



ei 



1 

n H — 
2 



+ e 2 



n + -) 3 + e 3 + 



e 4 



n + 



2/3 



+ . . . 



(14) 
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where the first few coefficients e, are given by 

= 8192000000000000 f * ^ 

61 48413843986233 m 2 V m 2 



223553096351603200000 (^J^) 2/3 5 1 / 3 /^\ 2 / 3 2 
62 ~ 2759589107215281 \ VJ 1 J 

vW ^ ^ 

22736856054709769417272276 m 2 , ™ 2 ^ 4 

e 3 = « -0.0466914 (17) 

486959893945285848925248675 fi fi K J 

9339788533820875 5^ f /ih 4 \ 1/3 



e 4 



131900122669056 23 y/7 371496873^ \ m 2 



+ 



8271957929254630511939617633398178757 /m 10 a/\ 1/3 



0.030669 ( ^ j + 0.00424238 ( ) . (18) 



999241702375726561994610281100000000 23 18574843653 \ 

+0.00424238 » 
rrr y \ // 5 ft 

Notice that altough the coefficients of higher order could be calculated as well, 
it is not necessary to do so, since they would only provide corrections that would be 
dominated by the error that we are making evaluating the integrals Ji. 

In Fig. El we plot the energy eigenvalues of the last column of Table III of [Zj, 
corresponding to ft — 1, m = 1/2, uj = 2 and \i = 8000 and compare them with the 
approximate results obtained by using eq. (JHJ). It is clear that our formula, although 
quite simple, works excellently. 

In Fig. El we display the error over the energy defined as E 



En F(e ~i) x 100 

as a function of the quantum number n. The boxes have been obtained using our formula 
eq. ([HI) and assuming ft = 1, m — 1/2, u> — 2 and /i = 8000. In this case E^ xac£) are 
the energies of the anharmonic oscillator calculated with high precision in last column 
of Table III of [7]. The jump corresponding to n = 25 is due to the low precision of the 
last value of Table III of [Zj. The pluses and the triangles have been obtained using our 
formula eq. (fTlj) (pluses) and eq. (1.34) of jHj (triangles) and assuming ft = m = u = 1 
and fj, — 4. In this case E^ xact > are the energies of the anharmonic oscillator numerically 
calculated through a fortran code. We can easily appreciate that our formula provides 
an approximation which is several order of magnitude better than the one of eq. (1.34) 
of 0. 

As a further check on the soundness of our approximation we consider the quantum 
zeta function for the quartic anharmonic oscillator |9"1 HOj: 

OO 1 
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Figure 2. Comparison between the numerical eigenvalues of and the results 
obtained using eq. 1)14(1. We assume h= 1, m= 1/2, u> = 2 and p = 8000. 



The exact value of Z(l) corresponding to setting uj = and (i — 4 is known to be 
32/3 /lx 

Z(l) = — -r 5 - ) « 3.63500364488 . (20) 

87T^ \3/ 

We have estimated Z(l) by using the first few numerical eigenvalues and then 
using for the remaining ones the analytical formula corresponding to the solution of 
eq. ()14j) . Using our formula eq. (|14jl only taking the first 4 eigenvalues numerically we 
have obtained the value 

Z 4 (l) = 3.635002 (21) 

which can be compared with the value 

Z uo (l) = 3.6350017 (22) 

obtained in ^U] by using a much larger set of numerical eigenvalues. 

We have also applied our method to the case of a anharmonic potential of the 
form V(x) = ux 2 /2 + px 6 /6. In this case we have used the WKB expansion to leading 
order and we have approximated with our method the integral J\ corresponding to this 
potential. In Fig. HJwe display the error H = — — ( ex l ct) — x 100 as a function of the 

energy. j^ xact ^ j s the exact value of the integral calculated numerically, whereas jf^ is 
the analytical approximant obtained with our method to order 10. Although the error 
increases with the energy, we notice that it is always smaller than 0.01%. As in the case 
of the quartic oscillator, our approximation works better in the region of low energy, i.e. 
opposite to the WKB method. 

By using the approximation for the integral J\ obtained with our method we have 
derived a very simple formula for the spectrum of the sextic oscillator to leading order 
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Figure 3. Error over the energy of the anharmonic oscillator with % = 1, m = 1/2, 
lj = 2 and /i = 8000 (set 1) and with % = m = uj = 1 and u = 4 (set 2). The boxes 
and the pluses have been obtained with our method, the triangles correspond to the 
error calculated using the analytical formula of 



in the WKB expansion: 



where 



E n = ~ v 3/2 ( 23 ) 

(ai (n+|)+a 2 -y/3i (n + ±) 2 + /3 2 (n + |) + ft ] 



ft 



20887573200705159168 66 v^ft pe ^P 5/6 

ati = g — « 18.46505979 - 

5057378805833647189 m2 ™ 5/2 

= 81296300711504421076V 

5057378805833647189 mu 2 moj 2 

4799197856362980042268993891174358974464 6^ h 2 p§ 

255770803856953672746878259 25727601721 m 5 cu s 

340.9584332 —t— 

m b U) iS 

339616486411617501798889917918008573952 6« ^lTfrpi 
255770803856953672 74687825925727601721 ml cj 6 

196.0168989- ' 



m 7 / 2 u; 6 

177921136833163170599693588290032640752 6_s pi 
25577080385695367274687825925727601721m 2 a; 4 

2/3 

12.64038572 P 



m 2 u; 4 



In Fig.EJwe compare the lower part of the spectrum obtained using eq. ()23j) with the 
exact eigenvalues, obtained numerically with a fortran code. Once again the agreement 
between the two is impressive. Fig. (jUJ) displays the error (in %) over the energy obtained 
using eq. 1)23)1 . The error is of the order of a few percents for low values of n, but 
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Figure 4. Error over the integral J\ for U = m = cj = p = 1. 
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Figure 5. Lower part of the spectrum of the sextic anharmonic oscillator. The 
circles are the exact eigenvalues calculated numerically and the pluses are the results 
of eq. We assume % = m = w = p = l. 



drops quite quickly, reaching a plateau of about 0.01% at large n. The presence of the 
plateau is a consequence of the error that we are making in the evaluation of J7i, as 
already discussed when we considered Fig. HJ Since the method that we are using is 
convergent ^ I2j , by working to higher orders one can arbitrarily lower the position of 
such plateau, to obtain the highly excited part of the energy spectrum with the desired 
accuracy. 

In conclusion we have shown that the method of QjJ |2] can be easily applied to 
calculate fully analytical expressions for the energy spectrum of quantum potentials 
within the WKB approximation. By applying our method to a sufficiently large order 
in 5 we are able to estimate analytically the WKB integrals with high precision, given 
the convergent nature of the expansion that we are performing. The application of this 
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Figure 6. Error over the energy obtained using eq. H23[) and assuming h = m = lo 
p=l. 



method to more general potentials in 1 and 3 dimensions is currently underway. 
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